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ABSTRACT 



We study the equilibrium of pressure truncated, filamentary molecular clouds that 
are threaded by rather general helical magnetic fields. We first derive a new form of 
the virial equation appropriate for magnetized filamentary clouds, which includes the 
effects of non-thermal motions and the turbulent pressure of the surrounding ISM. 
When compared with the data, we find that many filamentary clouds have a mass per 
unit length that is significantly reduced by the effects of external pressure, and that 
toroidal fields play a significant role in squeezing such clouds. 

We also develop exact numerical MHD models of filamentary molecular clouds 
with more general helical field configurations than have previously been considered. We 
examine the effects of the equation of state by comparing "isothermal" filaments, with 
constant total (thermal plus turbulent) velocity dispersion, with equilibria constructed 
using a logatropic equation of state. 

Our theoretical models involve 3 parameters; two to describe the mass loading of 
the toroidal and poloidal fields, and a third that describes the radial concentration 
of the filament. We perform a Monte Carlo exploration of our parameter space to 
determine which choices of parameters result in models that agree with the available 
observational constraints. We find that both equations of state result in equilibria that 
agree with the observational results. Moreover, we find that models with helical fields 
have more realistic density profiles than either unmagnetized models or those with 
purely poloidal fields; we find that most isothermal models have density distributions 
that fall off as r -18 to r~ 2 , while logatropes have density profiles that range from 

to r -18 . We find that purely poloidal fields produce filaments with steep radial 
density gradients that are not allowed by the observations. 

Key words: ISM: magnetic fields - ISM: clouds - MHD 



1 INTRODUCTION 

Observations have revealed that most molecular clouds are 
filamentary structures that are supported by non-thermal, 
small-scale MHD motions of some kind, as well as large 
scale ordered magnetic fields (cf. Schleuning 1998). Never- 
theless, virtually all theoretical models assume spheroidal 
geometry. While spheroidal models are a reasonable geom- 
etry for molecular cloud cores, these cannot adequately de- 
scribe molecular clouds on larger scales. The goal of this 
paper is to fully develop a theory for filamentary molecu- 
lar clouds including the effects of ordered magnetic fields. 
It is our intent that this work should elevate filamentary 



clouds to the same level of understanding as that enjoyed by 
their spheroidal counterparts (cf. review McKee et al. 1993). 
This is an important step in star formation theory because 
filamentary molecular clouds ultimately provide the initial 
conditions for star formation. A clear understanding of the 
initial conditions is necessary if we are to understand the 
processes by which clouds produce their star-forming cores. 

Ostriker (1964) investigated the equilibrium of unmag- 
netized isothermal filaments; he found that the density 
varies as ~ r -4 in the outer regions. However, this solution is 
much too steep to account for the observed density profiles in 
molecular clouds. For example, Alves et al. (1998, hereafter 
A98) and Lada, Alves, and Lada (1998, hereafter LAL98) 
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use extinction measurements of background starlight in the 
near infra-red to find r~ 2 density profiles for the filamentary 
clouds L977 and 1C 5146. 

Most theoretical models for self-gravitating filaments 
have featured magnetic fields that are aligned with the ma- 
jor axis of the filaments. The pioneering work by Chan- 
drasekhar and Fermi (1953) was the first to analyse the sta- 
bility of magnetized incompressible filaments with longitu- 
dinal magnetic fields. Stodolkiewicz (1963) developed a class 
of isothermal models in which the ratio of the gas to mag- 
netic pressure (/3) is constant. The magnetic field in these 
models simply re-scales the Ostriker (1964) solution; thus, 
the steep r -4 density profile is preserved. The structure of 
the Ostriker solution is unchanged by the addition of a uni- 
form poloidal magnetic field. The stability of such models, 
including the effects of pressure truncation, has subsequently 
been determined by Nagasawa (1987). Gehman, Adams, & 
Watkins (1996) have considered the effects of a logatropic 
equation of state (EOS) on the equilibrium and stability 
of filamentary clouds threaded by a uniform poloidal field. 
Unfortunately, these models possess infinite mass per unit 
length as a result. 

Observations suggest that some molecular clouds may 
be wrapped by helical fields (Bally 1987; Heiles 1987). There 
is also some observational evidence for helical fields in HI fil- 
aments towards the Galactic high latitude clouds (Gomez de 
Castro, Pudritz, & Bastien 1997). In fact, helical fields repre- 
sent the most general magnetic field configuration allowed if 
cylindrical symmetry is assumed. A few authors have previ- 
ously modeled filamentary clouds with helical fields. These 
models are similar to the Stodolkiewicz (1963) solution in 
that the magnetic pressure is proportional to the gas pres- 
sure, so that the density becomes a re-scaling of the Ostriker 
(1964) solution. 

Our analysis replaces the assumption of constant (3 
with the assumption of constant flux to mass loading for 
the poloidal (eg. Mouschovias 1976, Spitzer 1978, Tomisaka, 
Ikeuchi, and Nakamura 1988) and toroidal fields. We show 
that the magnetic field in this case has non-trivial effects on 
the density distribution, and if fact results in much better 
agreement with the available data. We also explore the role 
of the EOS by constructing models using both an "isother- 
mal" EOS, where the total (thermal plus non-thermal) ve- 
locity dispersion is assumed constant, and the pure loga- 
trope of McLaughlin and Pudritz (1996, hereafter MP96). 
The effects of pressure truncation play an important role 
in our analysis. By including a realistic range of external 
pressures, appropriate for the ISM, we show that the mass 
per unit length of our models is significantly decreased from 
the untruncated value. We also derive a new formulation of 
the virial equation appropriate for the radial equilibrium of 
pressure truncated filamentary equilibria with helical fields. 
We use this equation to compare our models with real fila- 
mentary clouds and to establish strong constraints on their 
allowed magnetic configurations. 

How would helical fields arise? All that is required is 
to twist one end of a filament containing a poloidal field, 
with respect to the other end. Even if molecular filaments 
form with an initially axial magnetic field, a helical field 
is plausibly generated by any kind of shear motion (such 
as subsequent oblique shocks, torsional Alfven waves, etc.) 
that twists the field lines. 



It is not the purpose of this paper to examine how he- 
lical fields could be generated. The main point of this work 
is that, having recognized that most molecular clouds are 
undoubtedly filamentary, magnetized, and truncated by an 
external pressure, it is of considerable importance to inves- 
tigate equilibrium models of molecular clouds that contain 
quite general helical fields and pressure truncation. We em- 
ploy two main approaches in our theoretical analysis. Firstly, 
we derive a general virial equation appropriate for pressure- 
truncated filamentary molecular clouds, which we use to un- 
derstand the roles of gravity, pressure, and the magnetic field 
in the overall quasi-equilibrium of filamentary clouds. Sec- 
ondly, we develop numerical MHD equilibrium models that 
can be compared with the internal structure of real clouds. 

Our virial analysis demonstrates that poloidal fields 
always help to support the gas against self-gravity, while 
toroidal fields squeeze the gas by the "hoop stress" of their 
curved field lines. Helical fields may either support or help 
to confine the gas, depending on whether the poloidal or 
toroidal field component is dominant. We show, in fact, that 
it is very difficult to understand observed clouds without the 
notion of helical fields and the confining hoop stresses that 
they exert upon their molecular gas. 

Having found evidence for helical fields from our virial 
analysis, we construct numerical MHD models of filamentary 
clouds in order to investigate the internal structure of mod- 
els that are allowed by the data. It is noteworthy that our 
isothermal models with helical magnetic fields always pro- 
duce density profiles that fall off as r -1,8 to r -2 , in excellent 
agreement with the data. We show that the toroidal field 
component is responsible for the more realistic behaviour, 
and that purely poloidal fields result in density profiles that 
fall even more rapidly than r~ 4 . We also consider the pure 
logatrope of McLaughlin and Pudritz (1996) as a possible 
effective EOS for the gas. We find that our logatropic mod- 
els have somewhat more shallow density profiles, but many 
are also in good agreement with the existing data. 

A brief outline of our paper is as follows. We first 
present the results of virial analysis of self-gravitating, pres- 
sure truncated, filamentary clouds containing both poloidal 
and toroidal field (Section g). In Section ||, we follow this up 
with a detailed analysis of the equations of magnetohydro- 
static equilibrium describing self-gravitating filaments and 
discuss important analytic solutions to these. A full numer- 
ical treatment of the equations is given in Section ^ where 
we also constrain our 3-parameter models with a wide va- 
riety of filamentary cloud data. We discuss these results in 
Section \3l and summarize in Section pi 



2 VIRIAL ANALYSIS FOR FILAMENTARY 
MOLECULAR CLOUDS 

In Appendix A, we use the tensor virial theorem to construct 
a scalar form of the virial theorem appropriate for pressure 
truncated filamentary clouds containining arbitrary helical 
fields. After carrying out the manipulations therein, we ob- 
tain 

= 2 J PdV-2P s V + W + M, (1) 
where the gravitational energy per unit length is given by 
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f 9$ 
W = - / P r -Q^dV. 



(2) 



and M is the sum of all magnetic terms (including surface 
terms) : 



M = h IbUv 



B 2 zS + B ls 

47T 



V. 



(3) 



This equation is appropriate for a non-rotating, self- 
gravitating, filamentary molecular cloud whose length 
greatly exceeds its radius. For the remainder of this paper, 
all quantities written with a subscript S are to be evaluated 
at the surface of the filament; thus we write that our fila- 
ment is truncated by an external pressure Ps at radius Rs- 
We further reserve calligraphic symbols for quantities evalu- 
ated per unit length; W is the gravitational energy per unit 
length since there are no external gravitational fields and 
V is actually the volume per unit length, or cross-sectional 
area nRs 2 , of the filament. As we shall now show, W can be 
evaluated exactly for a filament of arbitrary internal struc- 
ture and equation of state. The mass per unit length m of 
the filament is obstained by simply integrating the density 
over the cross-sectional area: 



(4) 



m = 2-7T / rp(r)dr. 



Poisson's equation in cylindrical coordinates takes the form 
1 d ( d$\ 



r dr V dr ) 



AnGp. 



(•») 



By integrating, we find that the mass per unit length interior 
to radius r can be written as 



m{r) = ^r — 



(6) 



2G r dr r " 

Using this result in equation ^, the gravitational energy per 
unit length can be transformed into an integral over the 
mass per unit length: 

W = -2G / m'dm = -m 2 G. (7) 
Jo 

It is remarkable that the gravitational energy per unit length 
takes on the same value regardless of the equation of state, 
magnetic field, or internal structure of the cloud. The only 
requirements are those of virial equilibrium and cylindrical 
geometry. McCrea (1957) gave an approximate formula for 
the gravitational energy per unit length as W = -am 2 G 
(where a is a constant of order unity) based on dimensional 
considerations; thus, our exact result gives a = 1 for all 
cylindrical mass distributions. 

By considering a long filament of finite mass M and 
length L, we find that the gravitational energy scales quite 
differently for filaments and spheroids: 



sphere 



GM 2 



3 GM 2 



(8) 



where a depends on the detailed shape and internal struc- 
ture for spheroids. It is of fundamental importance that the 
gravitational energy scales with radius for spheroids, but not 
for filaments. McCrea (1957) used this point to argue that 
filaments possess stability properties quite contrary to those 



of spheroidal equilibria. For spheroids, which best describe 
molecular cloud cores, the gravitational energy scales as 
~ R . As long as the core is magnetically subcritical, there 
always exists a critical external pressure beyond which the 
gravitational energy must dominate over the pressure sup- 
port. The equilibrium is unstable to gravitational collapse 
past this critical external pressure. On the other hand, the 
gravitational energy of a filament is unaffected by a change 
in radius. Thus, the gravitational energy remains constant 
during any radial contraction caused by increased external 
pressure. If the filament is initially in equilibrium, gravity 
can never be made to dominate by squeezing the filament; 
all hydrodynamic filaments initially in equilibrium are stable 
in the sense of Bonnor (1956) and Ebert (1955). 

In Appendix B, we consider the Bonnor-Ebert stability 
of magnetized filaments. Beginning with a discussion of uni- 
form filaments, we show that a uniform filamentary cloud 
with a helical field, that is initially in a state of equilibrium, 
cannot be made to collapse radially by increasing the ex- 
ternal pressure. We also give a more general proof which 
extends the argument to non-uniform filaments of arbitrary 
EOS. Thus, we conclude that all filamentary clouds, that are 
initially in a state of equilibrium, are stable against radial 
perturbations. 

The virial theorem for filaments (equation ^J) is best 
used to study the global properties of filamentary molecular 
clouds. It is useful to define the average density, pressure, 
and magnetic pressure within the cloud as 



(p) 



(P) = 



V 

V 



{Pmaa) - StvV 



B 2 z dV. 



(9) 



Quite generally, we may write the effective pressure inside 
a molecular cloud as P = a 2 p, where a is the total veloc- 
ity dispersion. We emphasize that all of our models take 
a to represent the total velocity dispersion, including both 
thermal and non-thermal components. It is particularly im- 
portant to note that when we describe an equation of state 
as "isothermal" , we really mean that the total velocity dis- 
persion is constant. The average squared velocity dispersion 
is defined simply as 



, 2, <P> _ J v ° 2 P d V 



(10) 



where the average has been weighted by the mass as in 
MP96. 

With the above definitions, we easily derive a useful 
form of our virial equation (equation ^) : 



Ps 
(P) 



1 - 



M 

\W\ 



(11) 



where M and W are the total magnetic and kinetic energies 
per unit length defined in equations || and ^| and m V ir is the 
virial mass per unit length defined by 



rriyij- — 



G ' 



(12) 



We note that tTT/yif IS analogous to the the virial mass 
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5R(a 2 ) 
G 



(13) 



normally defined for spheroidal equilibria. Using the defini- 
tion of the average magnetic pressure given in equation ^| 
we may write the total magnetic energy as 



M =2{{P mag ) - P mag , s )V. 



(14) 



Using this result, along with the expression for the gravita- 
tional energy per unit length (equation (?]) in equation [n], we 
obtain another useful for for our virial equation after some 
algebraic manipulations: 



Ps_ 

(P) 



\Pmag) Pmag ,S 

(P) 



(15) 



where P m ag,s is the total magnetic pressure evaluated at the 
surface of the cloud: 



B 2 z3 + B? 



iag,S — 



05 



(16) 



Equation |l5| makes two important points. First of all, 
the poloidal component of the magnetic field contributes to 
the magnetic pressure support of the cloud through (P mao ) . 
Secondly, the toroidal field enters into equation n5l only as 
a surface term, through P m ag,s, which helps to confine the 
cloud by the "pinch effect" well known in plasma physics. 

All magnetic fields, whether poloidal, toroidal, or of a 
more complex geometry, are associated with currents that 
flow within molecular clouds and the surrounding ISM. For 
a filamentary cloud wrapped by a helical field, the toroidal 
field component implies the existence of a poloidal current 
that flows along the filament. A natural question is whether 
a return current outside of the filamentary cloud completes 
the "circuit", or whether the poloidal current connects to 
larger scale structures in the ISM. The answer to this ques- 
tion will likely depend on the mechanisms by which fila- 
ments form, which might be addressed by future analysis. 
If the current returns as a thin current sheet flowing along 
the surface of the filament, the toroidal field at the surface 
would be nullifi ed, and so would its confining effects. As we 
show in Section 2.4, this would make the available data very 
difficult to understand, indeed. However, if the return cur- 
rent is diffuse and extended throughout the surrounding gas, 
as in the case of protostellar jets (Ouyed & Pudritz 1997), 
there would be a net magnetic confinement of the filament, 
which is consistent with the observations. 

In Appendix C, we derive the virial relations for filamen- 
tary molecular clouds analogous to the well known relations 
for spheroidal clouds (Chieze 1987; Elmegreen 1989; MP96). 
We show that the two geometries result in differences only 
in factors of order unity. Most importantly, we use our virial 
equation [i] to show that Larson's laws (1981) are also ex- 
pected for magnetized filamentary clouds of arbitrary EOS. 
The reader may consult Table fcll to compare expressions 
for m, R, (p), and E for spheroidal and filamentary clouds. 



2.1 Unmagnetized Filaments 

From equation jH], we see that unmagnetized clouds obey 
the following linear relation: 

Ps 



(P) 



= 1 



(17) 



This equation is exact for any unmagnetized filamentary 
cloud in virial equilibrium regardless of the underlying equa- 
tion of state or details of the internal structure. Since equa- 
tion [n] contains only quantities that are observable, we 
have derived an important diagnostic tool for determining 
whether or not filamentary clouds contain dynamically im- 
portant ordered magnetic fields. 

We can use equation ^ to obtain the critical mass 
per unit length for unmagnetized filamentary clouds. 
We consider a thought experiment in which mass is grad- 
ually added to a self-gravitating hydrostatic filament. As 
the mass per unit length increases, the compression due 
to self-gravity drives the filament to ever increasing inter- 
nal pressures, while the external pressure remains constant. 
This process continues until the cloud is so highly com- 
pressed that Ps/(P) 0, beyond which no physical so- 
lution to equation [l7| exists. By equation [n], this happens 
when m = m V i r ; thus, the virial mass per unit length plays 
the role of the critical mass per unit length for unmag- 
netized filamentary clouds. 

For a prescribed EOS, this procedure leads to an un- 
ambiguous determination of the value of m^. Depending on 
the EOS, the mass per unit length either approaches m,h 
asymptotically as Ps/{P) — * 0, or achieves mj, at some fi- 
nite radius where Ps vanishes. 

There is, however, one subtle point that needs to be 
made. For an isothermal equation of state, we can unam- 
biguously write that mj, = m V i r . However, the velocity dis- 
persion varies with density for non-isothermal equations of 
state. Thus, (a 2 ) and hence m V i r (by equation |l^) may vary 
as the cloud is compressed by self-gravity. The critical mass 
per unit length ran is the final value that m„; r takes before 
radial collapse ensues, while quantity that applies 

equally well to non-critical states. 



2.2 Magnetized Filaments 

When there is a magnetic field present in a filamentary 
molecular cloud, the critical mass per unit length m mag is 
significantly modified from the result obtained for unmag- 
netized clouds in the previous Section. Using the same ar- 
gument as presented above, a magnetized cloud achieves its 
critical configuration when Ps/{P) — > 0: 



M/W 



(18) 



where M and W are the total magnetic and gravitational 
energies per unit length given by equations [l4| and We 
recall that A4 may be either positive or negative, depend- 
ing on whether the poloidal or the toroidal field dominates 
the overall magnetic energy. In general, we find that poloidal 
fields increase the critical mass per unit length beyond m v i r 
for hydrostatic filaments, while toroidal fields reduce the crit- 
ical mass per unit length below m^. Physically, the reason 
for this behaviour is that the poloidal field helps to support 
the cloud radially against self-gravity, thus allowing greater 
masses per unit length to be supported. The opposite is true 
for the toroidal field component, which works with gravity 
in squeezing the cloud radially. 

We may constrain the critical mass per unit length for a 
filamentary cloud if we have additional information regard- 
ing the strengths of the poloidal and toroidal field compo- 
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nents. For a nearly isothermal EOS, the magnetic critical 
mass per unit length is given by 



rn h < 1 + 



{Pmag) 



mag ,S 



(19) 



Since molecular clouds are in approximate equipartition be- 
tween their magnetic and kinetic energies (Myers and Good- 
man 1988a,b, Bertoldi and McKee 1992), (P mag ) is not likely 
to greatly exceed (P). Therefore, it is unlikely that the mag- 
netic critical mass per unit length m mag would exceed the 
hydrostatic critical mass per unit length m./, by more than 
a factor of order unity. 



2.3 Surface Pressures on Molecular Filaments 

Molecular clouds are surrounded by the atomic gas of the 
interstellar medium (ISM). Like the molecular gas itself, 
the total pressure of the ISM is dominated by non-thermal 
motions. The external pressure is extremely important to 
our analysis since it both truncates molecular clouds at fi- 
nite radius and helps to confine the clouds against their 
own internal pressures (see equation [ll]) . Boulares and Cox 
(1990) have estimated the total pressure (with thermal plus 
turbulent contributions) of the interstellar medium to be 
on the order of 10 4 K cm . However, some molecular 
clouds are associated with HI complexes, whose pressures 
are typically an order of magnitude higher than the gen- 
eral ISM (Chromey, Elmegreen, & Elmegreen 1989). There- 
fore, we will be absolutely conservative by assuming that 
the external pressure on molecular clouds is in the range of 
10 4-5 K cm~ 3 . This assumption almost certainly brackets 
the real pressure exerted on molecular clouds by imposing 
the total (thermal plus turbulent) pressure of the ISM as a 
lower bound, and the pressure of large HI complexes as the 
upper bound. 

While the above pressure estimate is appropriate for 
most filamentary clouds, which are truncated directly by 
the pressure of the surrounding atomic gas, we note that a 
second type of filament exists, in which a dense molecular 
filament is deeply embedded in a molecular cloud of irregular 
or spheroidal geometry. The best example of this type of 
filament is the J -shaped filament in the Orion A cloud. In 
such cases, the external pressure must be estimated using the 
density and velocity dispersion of the surrounding molecular 



2.4 Comparison With Observations 

We have seen that the magnetic field affects the global 
properties of filaments only through the dimensionless virial 
parameter jM/|W|. The virial quantity A4/|W| provides a 
very convenient index of whether a cloud is poloidally or 
toroidally dominated and to what degree. For clouds with 
positive _M/|W|, the net effect of the magnetic field is to 
provide support and the field is poloidally dominated (cf. 
equation hj). WhenX/|W| is negative, the net effect of the 
field is confinement by the pinch of the toroidal field, and 
the field is toroidally dominated. Since M is directly com- 
pared to the gravitational energy |W|, the magnitude of our 
virial parameter provides an immediate indication of the im- 
portance of the ordered field to the dynamics of the cloud. 



In Figure [l], we have used equation |Tl| to draw contours of 
constant _M/|W as a function of m/m V i r and Ps/(P). The 
_M/|W| = (dotted) line represents all helical field configu- 
rations, including the unmagnetized special case, which have 
a neutral effect on the global structure of the cloud. Thus, 
we see that the diagram is divided into poloidally dominated 
(dashed lines) and toroidally dominated (solid lines) regions. 

Since both m/m V i r and Ps/{P) are observable quan- 
tities, we can constrain our models by locating individual 
filamentary clouds on this diagram. However, we must first 
compute m/m V ir and Ps/{P) for each cloud by the follow- 
ing steps. For each filament, we have found values for the 
mass, length, radius, and average linewidth from molecu- 
lar line observations in the literature (see table |l| for refer- 
ences). The mass per unit length m is obtained by dividing 
the mass of the filament by its length allowing for inclina- 
tion effects by conservatively assuming all filaments to be 
oriented within 45° of the plane of the sky. Since the emit- 
ting molecule (usually 12 CO or 13 CO) is always much more 
massive than the the average molecule in molecular gas, the 
observed linewidth must be corrected by applying Fuller and 
Myer's (1992) formula: 



Avl t = Av 2 oba +S\n(2)kT(= 

\m 



rn bs 



(20) 



where rn oos is the mass of the emitting species, m is the mean 
mass of molecular gas, and T is the kinetic temperature of 
the gas. For a normal helium abundance Y = 0.28, m = 2.33. 
We have assumed a temperature of 20 K for the gas; the 
exact temperature chosen makes only a small difference in 
Avtot since the turbulent component of the linewidth always 
dominates on any scale larger than a small core. The velocity 
dispersion may be obtained from equation EG by 



Avtot 
V8 In 2 ' 



(21) 



We identify a with {a 2 ) 1 ^ 2 as defined by equation 
since this velocity dispersion is obtained from an average 
linewidth for the entire cloud. With {a 2 } known, we compute 
m V i r from equation |l2], which directly gives us m/m V i r . Ob- 
taining the average radius Rs directly from maps, and hence 
the cross-sectional area V, the average density and internal 
pressure are then easily obtained using equations and ^| 
All that remains to deduce .M/[W| from equation O is to es- 
timate the external pressure. Most of the filamentary clouds 
in our sample are surrounded by atomic gas. Therefore, we 
conservatively assume the total external pressure to be in 



the range 10 4 - 10 b K cm , as discussed in Section In 



fact, the only exception is the J -shaped filament of Orion 
A, which is deeply embedded in molecular gas. In this case, 
we have estimated the external pressure from measurements 
of t he density and linewidth in the Orion A cloud (See table 
2.4 for references). 

Figure [j] demonstrates that most filamentary clouds re- 
side in a part of parameter space where 

0.11 ~ m/m V i r ~ 0.43 
0.012 ~ P S /{P}~ 0.75, 

(22) 

which is indicated by the shaded box in Figure [[]. Thus, 
we find that filamentary clouds range considerably in their 
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Cloud 




Region 


M 


L 

(pc) 


Rs 
(pc) 


a 

(kms -1 ) 


Rof. 


Notes 


LI 709 






140 


3 6 


0.23 


0.479 


o 


I 


L1755 






171 


6.3 


0.152 


0.526 


2 


1 


L171 


2- 


29 




219 


4.5 


0.156 


0.534 


2 


1 


DL 2 


1 




Taurus 


600 


6,1 


0.5 


1.08 


4 


4 


J-fil 






Orion 


5 x 10 3 


13 


0.25 


1.41 


1 


2,3 














0.35 


1.13 


5 




nf| 








1.55 x 10 4 


87.3 


2.25 


1.54 


3 


1 


SF |f 








3.65 x 10 4 


300 


2.25 


1.29 


3 


1 



a Dark lane in Taurus including B18. See Mizuno et al (1995) for a more detailed map. 

b The J'-shaped filament in Orion A. 

c Northern filament in Orion (See reference). 

d Southern filament in Orion (See reference). 

Table 1. We have compiled data on filamentary molecular clouds from several sources. 

References: 1. Bally (1987), 2. Loren (1989), 3. Maddalena (1986), 4. Murphy and Myers (1985), 5. Tatematsu et al. (1993) 
Notes 1. Little star formation. 2. Dense cores, star formation. 3. Deeply embedded in Orion A cloud. 4. Associated stars. 



Cloud 


m 

(Mqpc' 1 ) 


(Mqpc- 1 ) 


m/m vir 


(P) 

(10 4 if cm" 3 ) 


P.S 


Ps/(P) 


L1709 


35.9 


107 


0.34 


24.3 


3.2 


0.13 


L1755 


25.1 


129 


0.20 


46.8 


3.2 


0.068 


L1712-L29 


45 


132 


0.34 


82.4 


3.2 


0.038 


DL 2 


86.6 


547 


0.16 


63.6 


3.2 


0.050 


J-ffl. 


355 


925 


0.38 


1.77 x 10 3 


64.8 | 


0.037 




647 


590 


1.1 


1.05 x 10 3 


64.8 


0.062 


NF 


164 


1.11 x 10 3 


0.15 


12.1 


3.2 


0.26 


SF 


112 


777 


0.15 


5.8 


3.2 


0.55 



" Determined from Bally's (1987) density estimate and the 12 CO linewidth given by Maddalena (1986). 

Table 2. We have reduced the data of table [l] to obtain m/m v i r and P$ / (P) for each filament. We assume an external pressure P$ of 
-^q4.5±0.5 C m~ s for all filaments except the J'-shaped filament of Orion A, which is deeply embedded in molecular gas. We also assume 
that all filaments are oriented within 45° relative to the plane of the sky. We only give central values in the table, but the corresponding 
error bars are shown in figure hi 



virial parameters. However, it is remarkable that most of the 
clouds in our small data set appear to reside in the part of 
the diagram where M/\W\ < 0. Thus, our virial analysis 
infers that the magnetic field in at least several filamentary 
clouds is probably helical and toroidally dominated. Grav- 
ity and surface pressure alone appear to be insufficient to 
radially bind the clouds in our sample. While this means 
that filaments must be quite weakly bound by gravity, we 
note that similar results have also been obtained by Loren 
(1989b) and BM92. 

It is natural to wonder to what extent these conclu- 
sions could be affected by uncertainties in the observational 
results. The dominant sources of uncertainty in Figure [j] 
are probably the uncertainties in mass per unit length sur- 
face pressures. However, we have assumed very conservative 
ranges for the surface pressures and inclination angles of the 
clouds. Therefore, we do not believe that observational un- 
certainties can account for the helical fields that are required 
by our virial analysis. We also note that a more detailed 
model including rotation of the filament would necessarily 



lead to the same conclusion of a helical field. Since rota- 
tion would tend to support the cloud against gravity, even 
stronger toroidal fields would be required to confine the gas. 



3 EXACT MHD MODELS OF FILAMENTARY 
STRUCTURE 

The virial treatment of the previous section is perhaps 
the simplest and most illuminating way to understand 
the physics and global properties of filamentary molecular 
clouds. While the virial equations |l| and |l| are convenient 
to use, and are in fact exact expressions of magnetohydro- 
static equilibrium, the analysis can say nothing of the inter- 
nal structure of the clouds. This is the advantage of the exact 
analytic and numerical models developed in this section. 



© 1998 RAS, MNRAS 000, |-|l| 



Helical Fields and Filamentary Molecular Clouds 7 




m/m 



Figure 1. Helical field models are compared with the observed 
properties of real filaments. Curves are shown for various values 
of the virial parameter .M/[W|. Positive values, corresponding 
to the dashed curves, indicate that the poloidal field is dominant, 
while negative values, corresponding to solid curves, indicate that 
the toroidal field is dominant. The dotted line represents all solu- 
tions that are neutrally affected by the helical field (including the 
unmagnetized solution). The J'-shaped filament appears twice, 
because we have used two independent data sets in our analysis. 



3.1 The Poloidal and Toroidal Flux to Mass 
Ratios 

We postulate that the magnetic field structure corresponds 
to that of constant poloidal and toroidal flux to mass ratios 
T z and F^. The meanings of the flux to mass ratios are 
illustrated in Figure ^, and are defined in the following way. 
Consider a bundle of poloidal field lines passing through a 
small cross-sectional area of the filament 8V. The magnetic 
flux passing through the surface is B Z 5V, while the mass per 
unit length is pV. Thus, the ratio of the poloidal flux to the 
mass per unit length is 



P 



(23) 



Is there an analogous quantity for the toroidal component 
of the field? In fact the toroidal flux has been defined and is 
commonly used in plasma physics (Bateman 1978). Here, we 
consider a bundle of toroidal flux lines with cross-sectional 
area 5 A that form a closed ring of radius r centred on the 
axis of the filament. The mass enclosed by the ring is 2nr8A. 
Thus, we may define the toroidal flux to mass ratio (per 
radian) as 



i- oi = 



rp 



(24) 



The simplest field configuration is that of constant F z and 
r^,. We note that constant T$ results naturally if a filament 
of constant T z and length L is twisted uniformly through an 
angle <j>. Then 



Bcf, _ rip 
~B~ Z ~ T' 

which leads to the result 



(25) 



Poloidal flux/mass ratio 
(per unit length) 



Toroidal flux/mass ratio 
(per radian) 






mass: dM=2TtprdA 
flux: dO ( j ) = B ( j ) dA 



mass: dM=pLdA 
flux: dO z =B z dA z 

r = ^ 

Z p 



B 



L 



r p 



Figure 2. A schematic illustration of the poloida and toroidal 
flux to mass ratios introduced in equations and 





r<4 = 



(!) 



(26) 



We shall always assume constant T z and r</> for the remain- 
der of this paper. 



3.2 An Idealized Model: Uniform Magnetized 
Filaments 

As an iside, it is illustrative to consider uniform filamentary 
clouds are affected by helical fields of constant Y z and T^. 
In this simple model, the magnetic field is uniform within 
the filament, but drops to zero outside. By equation [m], the 
toroidal field B$ increases as ~ r within the filament and 



falls off as i 



in the external medium. Thus, the toroidal 



field is associated with a constant poloidal current density 
within the filament. 

With the assumption of constant density and the above 
definitions of Y z and r^, equation |l^ can be expanded to 



= 



2 Ps , 

(J h 

P 



rip 

8tt 



g n_ 

2 8tt 2 



(27) 



The critical mass per unit length is obtained by setting Ps 
0: 



2cr 2 + T'tp/^n 
G + rl/An* ■ 



(28) 



The effects of external pressure and the magnetic field 
are transparent in this simple model. Pressure and the 
poloidal field cooperate in supporting the cloud. On the 
other hand, the toroidal field enters into equation |2t] in con- 
cert with gravity. A filamentary cloud with a helical field 
would be confined jointly by gravity, external pressure, and 
the pinch of the toroidal field. Without prior knowledge of 
the field strength and direction (by molecular Zeeman and 
polarization observations), the cloud may appear to be un- 
bound by gravity alone. 
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3.3 General Equations for Magnetized 
Filamentary Molecular Clouds 

We consider the equilibrium structure of a non-rotating, 
self-gravitating molecular cloud with a helical field of con- 
stant flux to mass ratios F z and r^. We consider two pos- 
sible equations of state for the gas: 1) the "isothermal" 
equation of state P — a 2 p where a is the total velocity 
dispersion and 2) the "pure logatrope" of MP96 given by 
P/P c = 1 + Aln(p/p c ), where P c and p c are the central 
(along the filament axis) pressures and densities, and A is a 
constant. MP96 find A ~ 0.2 for molecular cloud cores. Al- 
though their analysis was based only on cloud core data, we 
shall assume that the same value of A might apply to fila- 
mentary clouds as well. We use these two equations of state 
because they probably bracket the true underlying equa- 
tion of state for molecular clouds; MHD cloud turbulence 
probably results in an EOS softer than isothermal (MP96; 
Gehman et al. 1996), while the pure logatrope is the softest 
EOS to appear in the literature. 

It is convenient to work in dimensionless units where 
density and pressure are scaled by their central values p c 
and P c - We further define the central velocity dispersion by 

a 2 = (29) 

pc 

A natural radial scale is then given by 



ld_ 

r dr 



' d * 
v dr 



and magnetohydrostatic equilibrium 
B 2 \ d , Id 



d_ 

dr 



P + 



8tt 



d $ 1 

^ dr r 2 dr 



r 2 B 2 
8tt 



(34) 



(35) 



In Appendix D, we construct the mathematical framework 
to solve these equations numerically for both isothermal and 
logatropic equations of state. We show that a solution to the 
dimensionless equations is characterized by three parame- 
ters, namely the flux to mass ratios T z and defined by 
equations |2^ and |24| and a third to specify the (dimension- 
less) radius of pressure truncation. We express this third pa- 
rameter as a concentration parameter defined by C defined 
as 



C = loe 



(36) 



where ro is the radial scale defined by equation [S3[ We note 
that our definition of C is analogous to the concentration pa- 
rameter C = log 10 (r t /Vo) defined for King models of glob- 
ular clusters (See Binney & Tremaine 1987). Our concen- 
tration parameter differs only in that the tidal radius r t is 
replaced by the pressure truncation radius, and our ro is 
smaller by a factor of 3. While we use C primarily as a the- 
oretical parameter, we note that it is in principle observable. 



r ° = T^-> ( 3 °) 
47rGp c 

which defines the effective core radius of the filament. Fi- 
nally, we may define natural scales for the mass per unit 
length and magnetic field: 



2 0"c 



Bo = P, 



1/2 



(31) 



Thus, all quantities are written in dimensionless form as 
follows: 



r 


= r/r 


P 


= P/Pc 


in 


= m/mo 


P 


= P/Pc 


a 


= <t/«Tc 


$ 


= <&K 2 


B z 


= Bz/Bo 




— B^/Bo 



(32) 

Hereafter, we will only ever refer to T z and F^, in their di- 
mensionless forms: 



(33) 



For brevity, we will drop the tildes for the remainder of 
this section and the next (except for where ambiguity would 
result); all quantities hereafter are understood to be written 
in dimensionless form unless otherwise stated. 

Our basic dimensionless equations are those of Poisson 




3.4 Analytic Solutions 

Before discussing numerical solutions, we derive a few spe- 
cial solutions that can be expressed in closed analytic form. 
Specifically, we discuss the unmagnetized isothermal solu- 
tion that was found by Ostriker (1964) (a brief derivation 
is given in Appendix D.). We note that that this solution 
is a special case of a more general magnetized solution ob- 
tained by Stodolkiewicz (1963); for brevity, we shall refer to 
this solution as the Ostriker solution for the remainder of 
this paper. We also find a singular solution for logatropic 
filaments. It is unlikely that either of these special solutions 
describe real filaments, which are probably magnetized and 
non-singular, but they do serve as important benchmark re- 
sults to compare with our more elaborate magnetized mod- 
els. 



3.4-1 The Ostriker Solution: Unmagnetized Isothermal 
Filaments 

The analytic solution for the special case of an unmagnetized 
isothermal filament is easily obtained using the mathemati- 
cal framework in Appendix D. It was first given by Ostriker 
(1964): 



Pc 



(l + rV8rg)»' 



(37) 



where we have restored the dimensional units. We note that 
the density decreases as ~ r~ 4 at large radii. That such steep 
density profiles have not been observed could be explained 
by three possibilities: 1) molecular clouds are not isothermal. 
A softer EOS would give a less steeply falling density at 
large radius; 2) real clouds contain dynamically important 
magnetic fields that modify the structure of the filament 
at large radius; 3) real filaments are always truncated by 
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external pressure. If the filament is truncated before the ~ 
r~ 4 envelope is reached, such steep behaviour would not be 



observed. We demonstrate in Section 4.1.2 that either of 
possibilities 1) or 2) can explain the observed properties of 
molecular clouds. 



3.4-2 Singular Logatropic Filaments 

Although we have been unable to find the analogue of the 
Ostriker solution for the logatropic EOS, we have been suc- 
cessful in finding a singular solution. For this model, we 
reinterpret p c as the density at some fiducial radius. We 
postulate a power law solution of the form 

p oc r a , (38) 

and find that a solution can only be obtained if a — — 1. 
The final solution with dimensional units restored is 



ra 



(39) 



It is useful to compare our solution with the singular logat- 
ropic sphere found by MP96: 



Psphe 



'2A 



(40) 



where we have rewritten their solution using our definition 
of ro- (Their definition of ro differs from ours by a factor 
of 3. Our definition is the customary choice for filaments.). 
It is remarkable that both singular logatropic spheres and 
filaments obey precisely the same power law. 



4 NUMERICAL SOLUTIONS 

We now t urn our attention to numerical solutions of equa- 



tion D20 using various values of the flux to mass ratios Y z 
and r^. Many of the solutions are shown out to very large 
radius but may be truncated to reproduce any desired value 
of the concentration parameter C defined in equation [}6| 

4.1 Numerical Results 



We have shown in Section 2.4 that many filamentary clouds 
are probably wrapped by helical magnetic fields. However, 
before con sideri ng the most general case of helical fields 
in Section 4.1.2, we first separately consider the effec ts of 
poloidal and toroidal magnetic fields in Section E^Lj . We 
note that purely poloidal fields are not allowed by our virial 
analysis, and purely toroidal fields are probably unrealistic. 
Neverthess, this is the best way to understand the roles of 
each field component in our more general helical field mod- 
els. 



Equation D2C gives the set of differential equations that 
we integrate to produce our models. The integration was 
done in a straightforward manner, using a standard Runge- 
Kutta method. 



4-1.1 Models With Purely Poloidal and Toroidal Fields 

Figures [| and ^ show the density and pressure profiles, the 
magnetic structure, and the mass per unit length for isother- 
mal and logatropic filaments threaded by a purely poloidal 



field. We have also included the velocity dispersion and av- 
erage velocity dispersion (given by equation ^) for the lo- 
gatropic equation of state. 

On each set of figures we have shown the density, pres- 
sure, and velocity dispersion structure of the unmagnetized 
solutions with dashed lines. For isothermal solutions, we 
have also drawn a line representing the asymptotic r~ 4 be- 
haviour of the Ostriker solution. Similarly, the r _1 singular 
solution has been included on density profiles for logatropic 
filaments. These power laws are meant as a guide in inter- 
preting the asymptotic behavior of the solutions; we find 
that they are obeyed at large radius for all unmagnetized 
filaments. 

Comparing the unmagnetized solutions (dashed lines) 
in Figures ^ and ^j, we observe that the density profile of the 
unmagnetized logatropic filament is slightly more centrally 
concentrated than the isothermal Ostriker solution, but falls 
off much less steeply at large radius. These figures also show 
that isothermal and logatropic filaments both tend to finite 
mass per unit length, although they differ in that isother- 
mal filaments approach the critical mass per unit length only 
asymptotically as their radii tend to infinity. As discussed 
in Section this limit represents the critical mass per unit 
length m,h (see Section ^7l]) beyond which no equilibrium is 
possible. For the isothermal filament, it is easy to show ana- 
lytically (from equation ^) and we verify numerically that 
nih = 87rmo, where mo is the mass scale defined by equa- 
tion ^l[ For the unmagnetized logatropic filament, we find 
numerically that ran = 185. 8mo. Logatopic filaments can 
support a greater mass per unit length for equivalent central 
velocity dispersion a c . This is easily understood since the av- 
erage velocity dispersion (a 2 ) 1 ^ 2 always exceeds the central 
value offering more turbulent support to the filament. 

Perhaps the most notable feature of logatropic fila- 
ments is that they "self-truncate" at finite radius and den- 
sity p = p c exp(— 1/A) where the velocity dispersion and 
pressure vanish. The logatropic EOS is designed to have a 
nearly isothermal core and a rising velocity dispersion out- 
side of the core radius. At some point, however, the velocity 
dispersion turns over and falls to zero. The velocity disper- 
sion could of course never vanish in a real cloud since all real 
clouds are truncated by finite external pressure. Whether the 
region of outwardly falling velocity dispersion actually falls 
within the pressure tru nca tion radius in fits to real clouds 
is addressed in Section 4.2, where we attempt to constrain 



our models using the observational results of table 2.4. 

Because of the way in which we have defined T z (equa- 
tion ^3j), B z is exactly proportional to the density. The 
toroidal field, however, shows a more interesting structure; 

always vanishes along the axis of the filament, as it must 
for the field to be continuous across the axis. We note that 
the logarithmic radial scale of Figure ^] makes the vanishing 
of at the axis difficult to see in some cases. It also is 
found to decay at large radius for isothermal filaments since 
B<f, — r^rp oc r -3 . Hence, there is a single maximum in 
the toroidal field structure. For logatropic filaments, the r _1 
asymptotic behaviour of the density implies that B<f, = T^rp 
tends to a constant value at large radius. Thus, the toroidal 
field in logatropic filaments lacks the local maximum found 
for isothermal filaments. 

The effects of poloidal and toroidal magnetic fields on 
the density structure are apparent in Figures |^ and ^ For 
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either EOS, the poloidal magnetic field supports the cloud 
and causes it to be more extended radially than the corre- 
sponding unmagnetized filament. Toroidal magnetic fields, 
on the other hand, pinch the filament to smaller radial ex- 
tent. Thus, we find our numerical results to be in agreement 
with the results of our virial analysis in Section []. 

It is significant that purely poloidal fields always steepen 
the outer density profile, while toroidal fields make it more 
shallow. For isothermal filaments, poloidal fields always re- 
sult in density profiles that are steeper than the r -4 be- 
haviour of the Ostriker solution; this is true of even logat- 
ropic filaments when the field is of sufficient strength. These 
steep density profiles have never been observed, so purely 
poloidal fields do not match the data. Thus, it seems likely 
that the field must have a toroidal component if a realistic 
density profile is to be achieved. 

We find that the magnetic field has a dramatic effect 
on the critical mass per unit length m ma9 of the cloud. 
For either EOS, a poloidal magnetic field increases mmag, 
since the poloidal field acts to support the cloud against 
self-gravity. The toroidal magnetic field works with gravity, 
thus decreasing the maximum mass per unit length that can 
be supported. These conclusions are in agreement with our 
virial results from Section 3. 

From Figures ^| and ^ , we note that all isothermal fil- 
aments that are unmagnetized or contain a purely toroidal 
field tend to the same mass per unit length m^. This is easily 
explained by our virial equation jl5| since the toroidal field 
always tends to zero at large radius for isothermal filaments. 
The critical mass per unit length m mag clearly cannot be af- 
fected, since the toroidal field only enters the virial equation 
through its surface value. This is not the case for logatropic 
filaments because B^ tends to a constant value at large ra- 
dius. 



4.1.2 Helical Field Models 



In Section 2.4, we provided evidence based on our virial anal- 
ysis that filamentary clouds likely contain toroidally domi- 
nated helical magnetic fields. At this point, we shall take a 
further step by comparing our exact MHD models with the 
observed properties of filamentary molecular clouds. As we 
have noted in Section ^, a numerical solution is completely 
determined by the choice of three dimensionless parameters; 
the flux to mass ratios F z , F^, and the concentration param- 
eter C. 

Although Rs can be observed with little difficulty, ob- 
taining an accurate value for C is difficult because of the 
uncertainty in the core radius ro. According to equation po| 
the core radius depends on both the central density and 
velocity dispersion along the axis of the filament, both of 
which might be quite uncertain. We can, however, estimate 
a rough upper bound to C using the data of table |!]. We do 
not presently know whether the central (axial) velocity dis- 
persions of filamentary clouds are dominated by non-thermal 
motions, as the bulk of the cloud certainly is, or if the ve- 
locity dispersions are thermal, as they are in many low-mass 
cloud cores. Nevertheless, we do know that that a c must be 
at least the thermal value, which is 0.23 km s _1 , assum- 
ing a temperature of 15K. Central densities are probably 
less than about I0 4 cm -3 , which is typical of a core. There- 
fore, equation implies that ro is probably not less than 



~ 0.04 pc. In Table [I], we find that Rs ~ 0.5 pc for most (but 
not all) of the filaments in our sample. Therefore, equation 
implies that most filamentary clouds should have con- 
centration parameters that are less than approximately 1.1. 
This estimate should be treated with caution, considering 
the uncertainties and generalizations in our calculation. In 
particular, we note that larger filaments, such as the North- 
ern and Southern Filaments in the Orion region (See Table 
[l]) have radii that are many times larger than the value that 
we used in our calculation and may, therefore, have concen- 
tration parameters that exceed our upper bound. 

Three observable quantities shall be required to con- 
strain our theoretical models. We have previously (Section 
2.4) found the virial parameters Ps/{P) and m/m v i r to be 
useful in showing that toroidally dominated helical fields 
play an important role in the virial equilibrium of filamen- 
tary clouds. We use these parameters, as well as a third 
parameter specifying the ratio of average magnetic to ki- 
netic energy densities to constrain our models. Accordingly, 
we define a virial parameter 



X 



M 



(41) 



where M and K are the average magnetic and kinetic energy 
densities within the cloud defined by 



M 
K 



J v (Bl + Bl)dV 
\{P){° 2 ), 



(42) 



and V is the volume of the cloud (not to be confused with 
V). Myers and Goodman (1988a, b) have provided consider- 
able observational evidence that the average magnetic and 
kinetic energy densities are in approximate equipartition, 
with M m K to within a factor of order 2. Therefore, we 
impose the auxiliary constraint that 



X « 0[1) 



(43) 



for filamentary clouds with realistic magnetic fields. This 
equipartition of energy has been explained by attributing 
the non-thermal motions within molecular clouds to inter- 
nally generated Alfvenic turbulence (BM92). Since super- 
Alfvenic turbulence is highly dissipative, the Alfven speed 
poses a natural limit for the non-thermal velocity disper- 
sion (BM92). Thus, we expect a w va for molecular clouds. 
Defining the average squared Alfven speed as 



J™v 2 A dm' 



it is easy to show that 



X 



3<<r») ' 



(44) 



(45) 



Therefore, X m 1 is a natural result for magnetized clouds 
supported against gravity by Alfvenic turbulence. In the 
analysis that follows, we assume that 



0.2 < X < 5 



(46) 



for all reasonable models and that 0.5 < X < 2 is appropri- 
ate for our most realistic models. 
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4-1.3 Monte Carlo Exploration of the Parameter Space 

In this Section, we perform a Monte Carlo sampling of our 
parameter space in order to determine which values of T z , 
r^,, and C result in models that obey all of our constraints. 
The Monte Carlo analysis is very straightforward. We simply 
assign random values to the three theoretical parameters 
and compute helical field models using the mathematical 
framework of Section ^5 . Once a solution has been obtained, 
we compute m/m v i r , Ps/{P), and X using equations [l^, [jj 
and from which we easily determine whether or not the 
solution obeys our constraints. Figure ^ shows the results 
of our Monte Carlo analysis for isothermal models, while 
Figure ^ shows the results for logatropic models. Each point 
in these figures represents a model that obeys our constraints 
on m/m V ir and Ps/(P)\ models that fall outside of these 
constraints have been discarded. 

The grayscale in Figures [|a and ^a represent different 
ranges for X. The most likely range of X, with 0.5 < X < 2 
is shown as the lightest coloured points. The next darkest 
gray dots represent a less likely, but still possibly allowed 
range, with 0.2 < X < 5, while the darkest gray dots repre- 
sent models that are outside of these ranges, and therefore 
have unrealistically large or small magnetic fields. It should 
be noted that there relatively little overlap between these re- 
gions; they map out quite distinct regions on the diagrams. 
From Figure ^a, we find that the allowed ranges for the flux 
to mass ratios are approximately 



5 ~ r d 



25 

8 



(47) 



for isothermal filaments with 0.5 < X < 2. A somewhat 
larger region of the parameter space is allowed for filaments 
with 0.2 < X < 5. Comparing with Figure we find that 
the allowed flux to mass ratios are somewhat larger for lo- 
gatropic filaments, where we find 



10 
if) 



r 
r z 



40 
20 



(48) 



when 0.5 < X < 2. We note that Y z is more tightly con- 
strained than F0 for both isothermal and logatropic fila- 
ments. 

In Figures ^b and ^|b, we plot the our magnetic param- 
eter X, for the allowed models, against the poloidal flux to 
mass ratio F z . We find that X has a very strong dependance 
on r z for both isothermal and logatropic models. Moreover, 
we find that there is no obvious correlation between X and 
r^,. Since we can regard X as nearly a function of F z alone, 
the auxiliary constraint on X directly constrains F z . It is for 
this reason that somewhat tighter contraints are obtained on 
r z in equations ^ and |l8| than on F^,. 

Figures and ^c show the dependence of the concen- 
tration parameter C on m/m V i r and Ps / (P) for models that 
are allowed by the observations. We find that C may range 
from to « 3 for isothermal models, but C ~ 1.7 for most 
solutions where 0.5 < X < 2. Moreover, we find that C 
correlates rather well with Ps/{P), with greater values of 
Ps I (P) corresponding to smaller values of C. We note that 
most filamentary clouds probab ly have C ~ 1.1, considering 
our discussion in Section 4.1.2. However, we do not enforce 



this upper bound as a rigid constraint, since further data on 
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Figure 5. Isothermal Models. We show the results of our Monte 
Carlo exploration for isothermal filaments. Each point on these 
figures represents a mode that obeys the observational con- 
straints given in equation 2^; thus, we determine which ranges 
of r z , r^, and C result in models that agree with the available 
observational data, a) (top) The grayscale represents different 
ranges for X, as defined in equation ttlj. The most realistic solu- 
tions, with 0.5 ~ X <~ 2 are shown as white dots. The medium 
gray dots have 0.2 ~ X ~ 5. The dark gray dots represent mod- 
els that are outside of these ranges with unrealistic magnetic field 
strengths, b) (middle) We show that X is determined mainly by 
r z . c) (bottom) We show the allowed ranges of the concentra- 
tion parameter C. The shading is the same as in a). However, we 
note that the darkest gray dots are mostly hidden "behind" the 
medium gray, in this case. 



14 Jason D. Fiege and Ralph E. Pudritz 




s lf ' •« > ? «' 




Figure 6. Logatropic Models. We show the results of our Monte 
Carlo analysis for logatropic filaments. The graphs are as de- 
scribed in the caption of Figure EL 



the central densities and velocity dispersions of filamentary 
clouds needs to be obtained in order to make our argument 
definitive. We find that C ~ 1 whenever P S /{P) ~ 0.25; 
therefore, isothermal filaments with C ~ 1 must be subject 
to external pressures that are at least one fourth of the mean 
internal pressure. Such relatively high external pressures are 
well within the range of pressures allowed by equation |22[ 
The concentration parameter C is much more restricted for 
logatropic models, where C may range only from aprroxi- 
mately 0.4 to 1.2. As a general trend, we find that C in- 
creases slightly with m/m V i r , and also as Ps/{P) decreases. 
This is a natural result, since filaments become more radi- 
ally extended, with greater C, as they become closer to their 
critical configurations with vanishing Ps/{P) and maximum 
m/m V ir- 



4.2 "Best-Fitting" Models For Magnetized 
Filamentary Clouds 

In Figure ^, we show 50 isothermal helical field models that 
span the range of parameters allowed by equations ^2] and 
. We see that our allowed models possess a number of very 
robust characteristics. Most importantly, we find that most 
of our isothermal models have outer density profiles that fall 
off as ~ r~ 1,8 to r~ 2 , with some of most truncated models 
having somewhat more shallow profiles. This is most clearly 
shown in Figure ^b, where we have plotted the power law 
index a = dlnp/dlnr as a function of the dimensionless ra- 
dius r/rrj. We observe that a becomes more negative with 
increasing radius, but that none of our models ever have den- 
sity profiles that are steeper than r~ 2 . Thus, we find that 
our isothermal helical field models have density distributions 
that are much more shallow than the r -4 Ostriker solution. 
This radical departure from the Ostriker solution is clearly 
due to the dominance of the magnetic field over gravity in 
the outer regions. The overall effect of the helical field is 
to modify the density structure of the Ostriker solution so 
that a much more realistic form is obtained. In particular, 
we note that A98 and LAL98 have recently used extinction 
measurements of background starlight in the near infra-red 
to show that two filamentary clouds, namely L977 and IC 
5146, have r~ 2 density distributions. Our helical field models 
have density profiles that are essentially the same as those 
obtai ned for models with purely toroidal fields in Section 
|4.1.l| . Therefore, we conclude that the outer density distri- 
bution is shaped primarily by the toroidal component of the 
field. We note, however, that the toroidal field is in fact 
much weaker than the poloidal field throughout most of a 
filamentary cloud. In all cases, the basic magnetic structure 
is that of a poloidally dominated core region surrounded by 
a toroidally dominated envelope, where the field is relatively 
weak. 

In Figure |^, we show a sample of 50 logatropic mod- 
els that are allowed by our constraints. The main difference 
between the logatropic models and the isothermal models 
shown in Figure ^ is that there is a much greater variety 
of allowed density distributions for the logatropes. We find 
logatropic filaments with density profiles as shallow as r _1 
and as steep as r~ 1 ' 8 . Unlike the isothermal solutions, a 
does not decrease monatonically. Rather, it usually reaches 
a minimum value somewhat less than -1 when r/ro ~ 1 to 3, 
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and increases at larger radii. The result is that the density 
distribution usually contains a small region where the den- 
sity falls quite rapidly, which is surrounded by an envelope 
with a more gentle power law. Many logatropic models have 
density profiles that are too shallow to explain the A98 and 
LAL98 data. However, we have also found many logatropic 
models that approach the observed r~ 2 profiles. The main 
difference between isothermal and logatropic models is that 
isothermal filaments produce a nearly "universal" r~ 1,8 to 
r~ 2 density profile, while logatropic filaments show a much 
larger range of behaviour. 



5 DISCUSSION 

We show, in Section that most of the filamentary molec- 
ular clouds in our sample have velocity dispersions that 
are too high for them to be bound by gravity and surface 
pressure alone. Thus, we find evidence that many filamen- 
tary clouds are probably wrapped by helical magnetic fields 
whose toroidal components help to confine the gas by the 
hoop stress of the curved field lines. It is important to re- 
alize that this conclusion is based on our virial analysis of 
filamentary clouds (Section ^) and is, therefore, independent 
of either the EOS of the gas or the mass loading of the mag- 
netic field lines. We find that all of the filamentary clouds in 
our sample have masses per unit length that are much lower 
than the critical mass per unit length for purely hydrostatic 
filaments, and that filaments are quite far from their critical 
configurations, where Pg/{P) — > (See equation ^.). Thus, 
star formation in filamentary clouds must involve fragmen- 
tation into periodic cores (cf. Chandrasekhar & Fermi 1953) , 
rather than overall radial collapse. We refer to the second 
paper in our series for a full analysis of this process. 

We construct numerical MHD models of filame ntary 
clouds, which we explore thoroughly in Section 4.1.3 using 



Monte Carlo techniques. We show that both isothermal and 
logtropic models are consistent with the available observa- 
tional data. We find that helical fields have profound effects 
on the outer density profiles of filamentary clouds. While hy- 
drostatic fila ment s have a density profile that falls off as r -4 
(see Section 3.4. 1 ) , nearly all of our isothermal models with 
helical fields have much more shallow density profiles that 
fall off as r~ 1,8 to r~ 2 . The toroidal component of the field 
is entirely responsible for these shalow density profiles. Lo- 
gatropic models show a a greater variety of behaviour, with 
density profiles ranging from r _1 to r -1 ' 8 . We note nearly 
all of our isothermal models and many of our logatropic 
models are in excellent agreement with the r~ 2 density pro- 
files of the filamentary clouds L977 and IC 5146, which have 
recently been observed by A98 and LAL98. Helical fields, 
rather than the nature of the EOS, seem to explain the data 
by controlling the density profiles. 

The confirmation of helical fields in filamentary clouds 
will ultimately require direct observations of the field struc- 
ture using sub-mm polarization and possibly molecular Zee- 
man observations. Unfortunately, most polarization obser- 
vations of molecular clouds have been carried out in the op- 
tical and near infra-red regions of the spectrum; Goodman 
et al. (1995) have demonstrated that such measurements 
are likely a poor indicator of magnetic field direction in cold 
dark clouds. More promising is the prospect of observing 




Figure 7. We show a sample of 50 isothermal models that span 
the range of allowed parameters given in equations ^ and ^fi| a) 
(top) We show the density profiles of the models. The dashed lines 
represent the r~ 4 density profile of the Ostriker solution and an 
r~ 2 profile, which is in agreement with the observed density pro- 
files of filamentary clouds (Alves et al. (1998), Lada, Alves, and 
Lada (1998)). b) (middle) We show the how the power law index 
d\n p/ d\nr behaves with radius, c) (bottom) We show the be- 
haviour of the poloidal and toroidal components of the magnetic 
field. 
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Figure 8. We show a sample of 50 logatropic models that span 
the range of allowed parameters given in equations and a) 
(top) We show the density profiles of the models. The dashed lines 
represent the r -4 density profile of the Ostriker solution and an 
r~ 2 profile, which is in agreement with the observed density pro- 
files of filamentary clouds (Alves et al. (1998), Lada, Alves, and 
Lada (1998)). b) (middle) We show the how the power law index 
dlnp/dlnr behaves with radius, c) (bottom) We show the be- 
haviour of the poloidal and toroidal components of the magnetic 
field. 



the thermal dust grain emission of dark clouds, which will 
hopefully become commonplace in the near future with in- 
struments like the SCUBA polarimeter. In emission, we can 
be assured that any polarization is due to warm dust grains 
within the cloud being observed. The observational verifi- 
cation of our helical field models for filamentary molecular 
clouds will need to rely heavily on such observations. 



5.1 Observational Signatures of Helical Magnetic 
Fields 

Helical magnetic fields present an interesting and unique 
polarization pattern. One might expect to find polarization 
vectors aligned with some average pitch angle of the mag- 
netic field. However, Carlqvist (1997) has modeled the polar- 
ization pattern of helical magnetic fields and demonstrated 
that this assumption is incorrect. The signature of a helical 
field, assuming that the polarization percent remains small 
and that the cloud is optically thin, is actually polariza- 
tion vectors that are either aligned with or perpendicular 
to the filament axis with a possible 90° change in orien- 
tation at some radius. The reason for this counterintuitive 
behaviour is that any line of sight through the filament in- 
tersects a given radius not once, but twice. Since the order 
of polarizing elements is unimportant in the limit of small 
polarization, the combination of any symmetric pair of po- 
larizing elements results in a cancelation of any oblique com- 
ponent of the net polarization. Although Carlqvist 's model 
was done for absorption polarimetry, the same reasoning 
holds for emission. The overall pattern will only reflect the 
dominant component of the field along any line of sight; 
thus, our models predict that the innermost (poloidally dom- 
inated) regions of filamentary clouds should be dominated 
in emission by polarization vectors aligned perpendicular to 
the filament, while the polarization should be predominantly 
parallel to the filament in the (toroidally dominated) outer 
regions. 

There is some direct evidence that filamentary clouds 
contain helical magnetic fields. Optical polarization and HI 
Zeeman data are consistent with a helical field in the large 
Orion A filament oriented at an approximately 20° pitch an- 
gle relative to the axis of the filament (Bally 1989) . While the 
reliability of the optical polarization data is uncertain due 
to the reasons discussed above, the Zeeman observations do 
not suffer such ambiguity. The line of sight component of the 
magnetic field has been observed to reverse across the L1641 
portion of the filament (Heiles 1987, 1989). Such a reversal is 
very suggestive of a helical field wrapped around the cloud. 
More recently, Heiles (1997) has suggested that an alternate 
explanation might be that the shock front of the Eridanus 
superbubble has swept past the Orion A cloud causing field 
lines to be stretched over the filament, thus simulating a he- 
lical pattern. While Heiles favours this idea, he cannot rule 
out the the older idea that the field is intrinsically helical. 
Regardless of the true nature of the J -shaped filament, we 
observe that it is one of the few filaments shown in Figure [l] 
where our models do not actually require a helical field for 
equilibrium. 
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6 SUMMARY 

1. All filamentary molecular clouds are truncated by exter- 
nal pressure that provides a total (thermal plus turbulent) 
pressure which is likely in the range of 10 4 to 10 5 Kcm~ 3 . 

2. We have derived a new form of the virial theo- 
rem appropriate for filamentary molecular clouds that are 
truncated by a realistic external pressure and contain or- 
dered magnetic fields. We have collected data on filamen- 
tary clouds from the literature. We find that most of the 
filamentary clouds in our sample are constrained by 



0.11 



~ m/m vir ~ 0.43 



0.012 ~ P S /(P}~0.75. 



(49) 

We use these observational constraints to show that many 
filamentary clouds are likely wrapped by helical magnetic 
fields. 

3. We have used our virial equation to derive virial re- 
lations for filaments that are analogous to the well-known 
relations for spheroidal equilibria (Chieze 1987; Elmegreen 
1989; MP96). We find that the virial relations for filaments 
differ from the corresponding relations for spheroids only by 
factors of order unity. 

4. We have studied the stability of filamentary molecu- 
lar clouds in the sense of Bonnor and Ebert. We find that 
all filamentary molecular clouds, that are initially in equilib- 
rium, are stable against radial perturbations. Thus, a cloud 
that is initially in equilibrium cannot be made to undergo 
radial collapse by increasing the external pressure; the only 
way to destabilize a filament against radial collapse is by 
increasing its mass per unit length beyond a critical value 
that depends on the magnetic field. This critical value is the 
maximum mass per unit length for which any equilibrium is 
possible; for purely hydrostatic filaments, the critical mas s 



per unit length is given by rrih = 2(cr )/G (See Section 2.1) 



The poloidal component of the magnetic field increases the 
critical mass per unit length by supporting the gas against 
self-gravity. The toroidal field works with gravity in com- 
pressing the filament; thus, the toroidal field decreases the 
critical mass per unit length. 

5. There are two exact analytic solutions that can be 
found in the limit of vanishing magnetic field. Stodolkiewicz 
(1963) and Ostriker (1964) have found the equilibrium so- 
lution for unmagnetized isothermal filaments; the density in 



this solution tends to an 



behaviour outside of the 



core radius. A98 and LAL98 have shown that the filamen- 
tary clouds L977 and IC 5146 have density profiles that fall 
off as r~ 2 . If this result holds true for other filaments as well, 
it is unlikely that the Ostriker solution describes real molec- 
ular filaments. In addition, we have found a singular solution 
for unmagnetized logatropic filaments. The density profile of 
this filament is much more shallow; it falls off only as r _1 , 
which is probably too shallow to agree with observational 
results. 

6. We have constructed exact nume rica l MHD models 
for filamentary clouds in Sections ^ and <L1. We have con- 
sidered both isothermal and logatropic equations of state, 
which likely bracket the true underlying EOS for filamen- 
tary molecular clouds. The magnetic field structure is more 
general than in previous studies; we have assumed only that 



the poloidal and toroidal flux to mass rat ios (T z and F^) are 
constant, which we justify in Section 3.2. 

7. Isothermal models with purely poloidal magnetic 
fields have density profiles that are even steeper than the Os- 
triker solution; it is unlikely that such models describe real 
molecular clouds. Toroidal fields result in density profiles 
that are more shallow than the Ostriker solution (typically 
r~ 2 profiles), and in better agreement with observations. 

8. We have performed a Monte Carlo analysis of our 
models, in which we randomly sample our parameter space 
and then determine whether or not the resulting model 
agrees with the observational constraints (equations and 
kq). We find both isothermal and logatropic filaments that 
are allowed by the data. We find that 



5 ~ r d 



25 

8 



for isothermal filaments, and 



10 ~ 
10 ~ r 2 



40 
20 



(50) 



(51) 



for logatropic filaments. 

9. Our best-fitting isothermal models have density pro- 
files that fall off as only ~ r - 18 to ~ r~ 2 , in contrast to 
the r -4 behaviour of the Ostriker solution. These shallow 
profiles are entirely due to the effects of the toroidal compo- 
nent of the magnetic field. Thus, helical magnetic fields are 
necessary for reasonable models of isothermal filaments. The 
logatropic filaments show a greater variet of density profiles 
that range from r~ x to r -1 ' 8 . Thus, some of the logatropic 
filaments may agree with the observed r~ 2 profiles (A98 and 
LAL98), while others may be somewhat too shallow. 
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APPENDIX A: DERIVATION OF THE VIRIAL 
EQUATION FOR FILAMENTARY 
MOLECULAR CLOUDS 

The tensor virial theorem is written in its most general form 
as 



1 d 2 I lk 

2 dt 2 



— 2T ik + Si, 



PdV + M v 



v 



+ 



47T 87T 



+ W lk - 2M lk 



Px k dS{M) 



(Chandrasekhar 1961) where we work in Cartesian coordi- 
nates Xi, and the integration is over the volume V of the 
cloud bounded by surface S. Here, P is the total pressure 
(with both thermal and non-thermal contributions) and I ik , 
T ik , Wik, and M ik are respectively the Cartesian tensors for 
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moment of inertia, kinetic energy, gravitational potential, 
and magnetic energy. The pressure and magnetic field both 
contribute bulk terms that support the cloud and surface 
terms that confine the gas. The magnetic term My is just 
the bulk component of the magnetic energy: 



M v = — I B 2 dV 

07T 



(A2) 



where B is the magnitude of the magnetic field. The tensor 
components are given by 



M ik 



dVpXiXk 



dVpViVk 

dVpx k —— 

axi 



871 



dVBiB k , 



(A3) 



where p is density, P is pressure, Vi and Bi are components 
of velocity and magnetic field, and $ is the gravitational 
potential. 

A filamentary molecular cloud may be idealized as a 
cylinder whose length greatly exceeds it radius. The most 
appropriate scalar virial equation for such an object is ob- 
tained by taking the sum of the x and y (1 1 an d 22) diagonal 



components of the tensor virial equation Al 



the 



z compo- 



nent can be ignored because we are concerned only with 
equilibrium in the radial direction. Considering the surface 
S to be the boundary of a filament of total volume V, we 
expand the ith diagonal component of equation Al using 



the definitions A3 



1 d?Iu 



2 dt 2 



PdV 



v 
1 

+ — 



PxidSi 



1 

47T 



+ 



1 

4tt 



(xiBiXBjdSj 



1 

87T 



9$ 
px t - — dV 

OXi 



B 2 Xi dSi, (A4) 



where summation over index j is implied, but no sum is 
to be taken over index i. We have also used the fact that 
the kinetic energy tensor Tik vanishes, since the velocity is 
everywhere zero. Summing over components 11 and 22 and 
setting In = I22 = for a filamentary cloud in equilibrium, 
we easily obtain 

= 2/ PdV — / Pr-dS- / pr-V^dV 
Jv Js Jv 



+ 



1-J^dV-^JjBUBDdV 



±.j s (r.B)(B. dS )-±.j s B>r.dS, 



(A5) 



where the position vector r and surface element dS are de- 
fined by 



r = xx + yy 
dS = dS x x + dS y y. 



(A6) 



Assuming that the magnetic field is helical, it is obvious that 
B ■ dS = 0; thus, the first magnetic surface term on line 3 



of equation CI is identically zero. Simplifying equation A5 



and dividing by the length, we obtain the form of the virial 
equation most appropriate for filamentary equilibria: 



(A7) 



(A8) 



(A9) 



= 2^ PdV-2P s V- J Pr^dV 

+ _L f Bldv- ( Bls + B * 3 \v. 

4-7T J V 4-7T 

It is interesting to note that this may be rewritten as 

Q=^u(l-^ \ +1V + .VI. 

where U is the internal energy per unit length: 

U = I J pa 2 dV. 

Thus, we see that our virial equation for filamentary equi- 
libria does in fact differ substantially from the usual form 
for spheroids. 



APPENDIX B: BONNOR-EBERT STABILITY 
OF MAGNETIZED FILAMENTS 

We have shown, in Section ^, that all hydrodynamic fila- 
ments, which are initially in equilibrium, are stable in the 
sense of Bonnor (1956) and Ebert (1955). In this section, we 
examine the Bonnor-Ebert stability of uniform magnetized 
filaments. 



Bl Uniform Clouds 

We begin by considering a radial perturbation of a uniform 
filament in which the mass per unit length and the velocity 
dispersion are conserved. Solving equation for the exter- 
nal pressure, we easily obtain 



Ps = 



— + — 

V V 



(Bl) 



where 



2 8f 



8tt 



(B2) 



Differentiating with respect to V and simplifying using equa- 
and ^| we obtain 



B2 



--(Ps + P mag ) <0 



(B3) 



tions P3 1 , 

dP s 
dV 

for all choices of V and Ps- We conclude that all uniform 
filaments in equilibrium (having m <= m mag ) are stable in 
the sense of Bonnor and Ebert. 



B2 Non-Uniform Clouds 

We now extend the argument to general magnetized filamen- 
tary equilibria using the mathematical framework of Section 



3.3 



and Appendix D. 
Again we consider a radial perturbartion of a filamen- 
tary cloud in which the mass per unit length m is conserved. 
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(B4) 



Following MP96, we shall also require that the central veloc- 
ity dispersion o c remain unchanged. Referring to equations 
|3l| and |3^, we may write 

m AirGm 
m = = = — . 

m at 

Thus, we see that the dimensionless mass per unit length 
m is also conserved during the perturbation. Since rh is im- 
plicitly a function of r alone, the dimensionless radius Rs 
must remain fixed during the perturbation. Therefore, we 
find that the radial perturbation takes the form of a sim- 
ple rescaling of the dimensionless solution with none of the 
dimensionless variables perturbed whatsoever. With this re- 
sult in hand, we write 



Ps = PcPs = otp c P s 



-1/2 



Rs, 



(B5) 



where a c , Ps, and Rs all remain fixed during the perturba- 
tion. Eliminating p c , we obtain the result 



— 1/2 

Since Rs oc Pa , we find that 



dRs 
dP s 



< 



(B6) 



(B7) 



for all external pressures Ps. Therefore, we conclude that 
all self-gravitating filaments that are initially in a state of 
equilibrium (which requires m < m mag by equation jl8|) , are 
stable in the sense of Bonner and Ebert. 



APPENDIX C: VIRIAL RELATIONS AND 
LARSON'S LAWS FOR FILAMENTS 

We derive the virial relations for filamentary molecular 
clouds analogous to the well known relations for spheroidal 
clouds (Chieze 1987; Elmegreen 1989; MP96). Using equa- 
tion [7] in our virial equation for filamentary clouds (equation 
0), we write 



= 2m(cr 2 } - 2P S V - am 2 G + M 



(CI) 



We have shown in Section |2| that the constant a is exactly 
unity for any filament but we retain the constant in order to 
more directly compare with the corresponding expressions 
for spheroidal clouds. Equation ^l] can be rewritten as 

1-M/\W\ 



TTlyi 



1-Ps/(P) m 



(C2) 



where we have introduced the observable virial parameter 
Q-mag for filamentary clouds analogous to that of Bertoldi 
and McKee (1992, hereafter BM92). We may also write 



M 

IWI 



(C3) 



where a non 
limit: 



is just ctmag evaluated in the unmagnetized 

(C4) 



1 - P S /{P) ■ 
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Tabl e C I. The coefficients of the virial relations given in equa- 
tions |C5| , We have included a correction sec i for the projected 
surface density of filamentary equilibria and defined the effective 
mass per unit length for a spere as m ap h ere = M/(2R). 



Equations C2 and C4 are easily solved for the mass per 
unit length m, radius Rs, average density (p), and surface 
density E. In order to compare with the corresponding re- 
sults for spheroidal magnetized equilibria (See MP96), we 
present the results as general expressions with the coeffi- 
cients written in table CI: 



R 

(P) 
E 



( in 

Cr 
C, 



V 2 1 (a 2 ) 

/2 



(GPs) 1 

x maq v ^' 



= C E - 



, / anon \ PS 

9 \ a n0 n — a) (cr 2 ) 

1 / anon A 1 / 2 (PS\ 1/2 

,ll 2 \anon-a) \G) 



(C5) 



We have included a correction for the inclination i of 
the filament the the plane of the sky; this only affects the 
expression for the surface density. These expressions are re- 
markably similar to those for spheroidal clouds, retaining 
identical functional forms and differing only by coefficients 
of order unity. Of course, the mass per unit length of a fila- 
ment cannot be compared directly to the mass of a spheroid. 
However, if we define an effective mass per unit length by 
taking m ap here = M/(2R), the resulting expression again 
differs from ours by only a numerical factor of order unity. 

As long as the external pressure is approximately con- 
stant, we find that a oc R 1 ^ 2 and E oc const x seci for 
filaments. Thus, Larson's laws apply to filamentary clouds, 
aside from a trivial geometric correction factor to take into 
account the inclination of the filament with respect to the 
observer. 



APPENDIX D: MATHEMATICAL 
FRAMEWORK 

In this Appendix, we construct the mathematical framework 
used to compute the numerical solutions described in Section 
^| We define the magnetic pressure of the poloidal field P mag 
and a useful quantity that depends on the toroidal field 
by 



mag 

b,i, 



r^/(87r). (Dl) 
Defining the effective enthalpy of the turbulent gas by 
dP 

dh=—, (D2) 
P 
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(D3) 



and introducing the magnetic potentials f z and / 
dP m ag 

df z = 

P 

dU = ~o~ . 

r 2 p 

we may write equation ^ as 

J;(fc + * + /„ + /#)=0. 

We are free to specify boundary conditions for each of these 
potentials along the filament axis; h = $ = f z — f$ — at 
r — 0. Thus, equation D4 can be integrated: 



(D4) 



h + $ + U = 0. 



(D5) 



It is useful to define a new radial variable by the trans- 
formation 



s = In (r/a), 



(D6) 



where a is a constant that is nec essary only to derive a 



special analytic solution in Section 



3.4.1 



Under this trans- 



formation, equations ^ and D4 can be combined to give 

£- 

ds 2 
where 



(h + fz + U) 



(D7) 



2 2 2s 

r p = a e p. 



(D8) 



To solve the differential equation D7, we must first ex- 



press h, f z , and f$ in terms of s and the new quantity <]/. 
Applying equation |d^ to the isothermal and logatropic equa- 
tions of state and integrating, we obtain the following for- 
mulas for the enthalpies: 



hiso = hip = In ^ — 2 In a — 2s. 
in,,,, = :il - J) .1(1- .>V"v|; : ) 



(D9) 



We have also used the definitions of s (equation D6) and \& 
(equation D8) to write the final forms of the enthalpies in 
terms of these variables. 

Assuming constant flux to mass ratios T z and T^, we 
derive from equations |5^, and 



Dl 



-p2 2 

8tt 



rjy 2 

8Tva 4 e 4a 
T 2 ^ 2 



8tt 



(D10) 



Substituting these relations into equations D3 and integrat- 
ing, we obtain 



4-7T 
£ 2 

47T 



(P-1) 



El 

4% 



(— 

V a 2 e 2 



(Dll) 



where we have applied our boundary conditions that both 
f z and fif, vanish along the axis of the filament where p = p c 
and * = 0. We have expressed all quantities in equations 
D9|; 



and 



Dll 



in terms of s and \l/ alone. Thus, equation Dl is 



closed and can, at least in principle, be solved for Since 



and are written in terms of $ (equations D8 



and D1C), these quantities may be determined. Finally, the 



poloidal a nd t oroidal magnetic fields may be obtained from 
equations 



Dl 



3.4.1 



As a brief example, we show how equation D7 naturally 
leads to the Qstriker (1964) solution discussed in Section 
From equation D7, we easily obtain 



& 2 ^ 



(D12) 



for isothermal equilibria. This equation can be solved in 
closed form; the solution is simply 

* = 2 sech 2 s. (D13) 
Converting back to r and p, the equation takes the form 

P = (l + //<*? ■ (D14) 

The boundary condition at r = is p — 1 in our dimen- 
sionless units; therefore, we require that a = y/S, giving the 
Ostriker solution of equation ml. 

It is most convenient for numerical solutions to write 



equation D7 as a pair of first order equations. This is best 



accomplished by writing the gravitational acceleration as 

(D15) 
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d a 1 d a 
-$ = — <E>. 



dr r ds 

Then Poisson's equation ^ becomes 
Id 

Note that there is no reason for numerical solutions to retain 



(D16) 



the constant scale factor a introduced in equation D6, For 



the remainder of this section, we will take a = 1. Wit h th e 
help of equation Dl and the definition of *t (equation D8), 
our numerical system becomes 



ds 



(h + fz + U 



rg 



d_ 

ds 



(D17) 



wher e h, f z , and f$ ar e fun ctions of ^ by equations D9 



and 

Dll. Thus, equations D17 can be rewritten as differential 
equations for \& and g. We write these equations explicitly 
below. 

Denoting t he d erivative 4- with a prime ('), we derive 
from equations D3 

F 2 

fz = —e 
1 4tt 

r 2 

* - 4^' 
Using equations 



3 (*' -2*) 



(D18) 



D9 



we must calculate h' separately for the 
isothermal and logatropic equations of state. We note that 
h! can be written in the general form 



ti = J ffi(s,*) + if 2 (s,*)*', 



(D19) 



where the functions H\ an d H 2 are given in table Dl 
and 



Using e quati ons D18 
of equations 



D19, we express our system 



in its final form: 
2s \]> 



Hz 



_± 1 i^ p -2s 



-(g + e- 3 *). 



(D20) 
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Hi 



Hz 



Isothermal 
Logatropic 



-2 



$-1 



Table Dl. Hi a nd are the functions related to the enthalpy 
by equation |D1S|. 



Since H2 is positive definite, this dynamical system is regular 
on the entire interval s 6 (— 00,00). 

These equations are now in a form that can be numer- 
ically integrated given appropriate initial conditions at the 
axis of the filament (r = 0). The problem arises however 
that r — occurs at s = —00 in our transformed variable; 
thus, we start the integration at a small but finite value of r. 
We expect that p tends to a constant value of unity near the 
axis for any no n-si ngular distribution. From the definition 
of (equation D8), we find that ^0 ~ e s ° where so is the 
initial value chosen for s (typically w —10). Recalling that 
g is the gravitational acceleration, we apply Gauss's law to 
find the initial value for g: 

So = ~ = -ie so . (D21) 
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